Condensed Matter Physics, 1998, Vol. 1, No 3(15), pp. 453-462 



< 

OO 
(N 






o 
en 



THIRD-ORDER RELATIVISTIC 
en: DYNAMICS: CLASSICAL SPINNING 

S: PARTICLE TRAVELLING IN A PLANE 

(N 

O-t- Roman Matsyuk 



Institute for Applied Problems in Mechanics and Mathematics, 

National Academy of Sciences. 

15 Dudayev Str., 290005 L'viv, Ukraine 

matsyuk<3l ms.lviv.ua 

Received March 12, 1998 



Mathisson's 'new mechanics' of a relativistic spinning particle is shown to fol- 
low, in the case of planar motion, from only general requirements of relativistic 
invariance and of the dependence on third order derivatives along with the 'vari- 
ationality' feature. The hamiltonian counterpart ultimately recovers the Dixon 
^ , equations for this case with the Pirani supplementary condition. 

Q^ ' Key words: Lagrangian, Hamiltonian, Ostroghrads'kyJ_mechanics, 

■"nI" . Classical^spin, Relativistic-top. 



PACS: 03.20. +i 02.40. +m 



1. INTRODUCTION 



^ \ Following the updating tendencies in the formal theory of variational calculus 

j^ \ promoted by the development of the intrinsic differential geometry as well as global 

analysis, many authors have revisited the Ostroghrads'kyj mechanics with higher 

derivatives. The subject develops continuously and, surprisingly enough, models of 

physical meaning breed (see [1] and references in [2]). 

In this paper we consider a model third-order dynamics of a classical parti- 
cle which, although restricted to the unrealistic three-dimensional flat space-time, 
provides an instructive example of how new hamiltonian systems of physical mean- 
ingfulness may arise from higher-order variational calculus. This example admits 
a comprehensive solution of the involved variational inverse problem for invariant 
third-order equation of motion. It turns out that the equation thus obtained may 
be interpreted as yet another description of a planar motion of the classical spin- 
ning particle in special relativity. Taking into account that the general-relativistic 
equation of motion of the gravitational dipole particle admits, among others, also a 
solution of only two degrees of freedom (see [3,4]), we hope that the results of present 
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investigation may contribute to the future lagrangian and hamiltonian formulation 
of the general Mathisson equation [5] in the realistic four- dimensional curved space- 
time: 

UqU'' = 1 . 

If only a free and planar motion is going to be considered, the equation (1) splits 
into the following two {fi, z/ = 0, 1, 2): 

moiif, - Sf.^u" = , (2) 

©3,.M^ = . (3) 

We shall demonstrate later that the equation (2) may be cast into the form of 
being the only third-order relativistic equation admitting a lagrangian description. 

Applying to it a kind of the 'hamiltonization' prescription of [6] then yields an 
equivalent to the Dixon equations [8] adopted to the case presently considered here, 
from which the equation (2) follows in turn, provided the Pirani supplementary 
condition 

Uq©P^ = (4) 

is in force. 



2. LAGRANGIAN DESCRIPTION 

Necessary and sufficient conditions to the existence of a Lagrange function for a 
third-order differential equation 

A . v" + (v'. ^v) A . v' + B . v' + c = (5) 

were established in [9]. They are expressed by means of the following system of 
partial differential equations in the independent variables t, x", and v" 

(i) djaf^bc] = 

(tt) 2B[,b]-3D^A„fe = 

(iii) 2 d^iaQb]c - 4 d^ia^b]c + d^c Kb + 2 D^ d^, Aab = 

(iv) d^(aCb) - D^ B(„6) = 

(V) 2 d. djaCb] - 4 djaQb]c + D 2 d Aab + 6 D 9 [^Afe,] = 

V V X. J_ <• I X. 

4 d^iaCb] - 2 D^ d^iaCb] - D^ -^ A,, = . 



(6) 



[VI 



454 



THIRD-ORDER CLASSICAL SPIN DYNAMICS IN PLANE 



We recall that the skew-symmetric matrix A, the matrix B, and the column 
vector c depend on the variables t, x, v = dx/dt. The differential operator D^ 
denotes the first order generator of the Cartan distribution, 

Any third-order Euler-Poisson equation (i.e. that of some variational origin) in two 
space dimensions fits into the form (5). 

We are interested only in the equations bearing the Poincare symmetry with the 
infinitesimal generator X parametrized by means of a skew-symmetric matrix f2 and 
some vector tt: 

X= — (tt • x) dt + goot 7z . dx + ^ • (y. a dx) 
+ 5fuu TT . 9v + (tt • v) V . 9v + r2 • (v A 9v) 
+ 2 (tt • v) v'. d^' + (tt • v') V . 9^' + ri • (v' A d^') . 

The centered dot symbol denotes the inner product of vectors or tensors; the 
lowered dot symbol denotes the contraction of a row-vector and the subsequent 
column-vector, or (sometimes) the contraction of a matrix and the subsequent 
column-vector. In order to chose a convenient expression of the symmetry concept 
we introduce a vector differential form e, associated with the equation (5): 

ea = A,5dv'^ + k,dt , 

k = (v' • dy) Av' + Bv' + c . (7) 

Now it is possible to cast the idea of the symmetry of the equation (5) into the 
framework of the concept of exterior differential system invariance. The system in 
case is generated by the vector valued Phaff form e and the contact vector valued 
differential forms 

dx - vdt , dv - v'dt . (8) 

Let X{e) denote the Lie derivative of the vector valued differential form e along the 
vector field X. The invariance condition consists in that there may be found some 
matrices ^, S, and 11 depending on v and v' such that 

X(e) = # . e + H . (dx - vdt) + n . (dv - v'dt). (9) 

We also assert that A and k in (7) do not depend neither on t nor on x. 

The identity (9) splits into more identities, obtained by evaluating the coefficients 
of the differentials dt, dx, dv, and dv' independently: 

(tt . ^v + (tt • v) V . ^v + ri • (v A ^v)) A 

+ 2(7r • v)A+ (Av) (8)7r - Ari = $A; (10) 

2 (Av') TT + (tt • v') A = n ; (11) 

-k(g)7r = H; (12) 

X(k) = #k - 3v - nv' . (13) 
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In the above the '0' symbol means the tensor (sometimes named as 'direct') product 
of matrices. 

A skew-symmetric two-by-two matrix always has the inverse, so the 'Lagrange 
multipliers' $, S, and 11 may explicitly be defined from the equations (10-12) and 
then substituted into (13). Subsequently, the equation (13) splits into the following 
identities by the powers of the variable v' and by the parameters O and tt (take 
notice of the derivative matrix A' = (v'. dy) A; also the vertical arrow sign points 
to the very last factor to which the aforegoing differential operator still applies): 

(fl • (v A ^v)) AV + {fl ■ (v' A ^v)) Av' - (v'. ^v) Afiv' 

= (O • (v A ^v)) AA-^AV - Af2A-^AV ; (14) 

(f2 • (v A ^v)) B - Bri = (f2 • (v A ^v)) AA-^B - AfiA'^B ; (15) 

(17 • (v A ^v)) c = (17 • (v A ^v)) AA-^c - AflA'^c ; (16) 

(tt . ^v + (tt • v) V . ^v) A'v' + (tt • v) AV + (tt • v') (v . ^v) Av' + (tt • v') A'v 
= (tt . ^v + (tt • v) V . ^v) AA-^AV + (ttA-^AV) Av - 3 (tt • v') Av' ; (17) 

(tt . ^v + (tt • V) V. ^v) B + (BV) 77 

= (tt. ^v + (tt • v) V . ^v) AA-^B + (Av) ® ttA'^B + (tt • v)B ; (18) 

(tt . 5v + (tt • v) V . ^v) c 

= (tt . ^v + (tt • v) V . ^v) AA-^c + 3 (tt • v) c + (7rA"^c) Av . (19) 

Cumbersome although routine calculations accompanying the simultaneous solv- 
ing of the partial differential equations (14) and (17) with respect to the unknown 
function A12 produce the unique output of 

const 



^12 



;i+ViVi+V2V2)3/2 • 

We remind that the system of the equations (14-19) and the system (6) must be 
solved simultaneously. Thus, the equation (6i) becomes trivial now. 

Under the assumption of B being a symmetric matrix (see (6ii)), the solution of 
the equations {(15), (18)} amounts to 

Bab = const ■ (1 + \l-v)-^^^[Vayb " (1 + V-v) gab] ■ 

This automatically satisfies the equation (6iii) too. In what concerns the subsys- 
tem {(16), (19)}, only the trivial solution c = exists. 

We are ready now to formulate the summary of the above development in terms 
of a proposition: 

Proposition 1 The invariant Euler-Poisson equation of a relativistic planar mo- 
tion is: 

■ " *^ 'v-v') + - ^ [(l+v-v)v'-(v'-v)v] = 0. (20) 



(l+Vv)3/2 (l+Vv)5/2^ ' (l+Vv)3/2 
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Arbitrary constant /i serves to parametrize the set of all the variational equations 
(20). The definition of the 'star operator' is common. Thus, *1 = e(i) A e(2) whereas 
* (e(i) A e(2)) = 1 if the (pseudo)orthonormal frame {e(i) , 6(2)} carries the positive 
orientation; also (*w)a = e^aW^. We recall for future use the definition of the inner 
product of two bi-vectors: 

(aAb)-(cAd) = (a-c) - (b-d) . 

Proposition 2 The Euler-Poisson equation (20) describes the free motion of a spin- 
ning particle in two space dimensions if 

u^S^"' = Q. (21) 

Demonstration. Equation (20) describes the world line of a particle parametrized 
by time. Passing to the proper time parametrization one obtains: 

u X u-\- ^ii = Q . (22) 

Let us introduce a vector a^ = ^e^XfiS'^^- Then the equation (2) takes the form 

m^,u -\- a X u = (23) 

with the consequence that 

a-u = 0. (24) 

Any vector a may always be presented as 

a ^ {a • u) u — {a X u) X u . (25) 

If the solutions of the equation (22) are to satisfy also the equation (23), we may 
define the variable u from (22) with the help of (it • ii) = — (■it • u) and substitute 
it into (23) to obtain, in view of (24), 

—fi (a • u) u — rrioU + (it • it) a x u = . (26) 

The condition (21) is equivalent to a x u = 0, thus the equation (26) gives 

(a- u) = —mo/fj,, (27) 

and from (25) it follows that 

a = -u, 28 

so the equations (22) and (23) are now equivalent. 

Remark 1. In most general setting, when ||u|| ^ 1 and |0| 7^ 1 (g = det(0^„), m, 
n run from to 3) the four-vector of spin is introduced by 

5p = 4^emnqpU"^©"V (29) 
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It is straightforward that S3 = —a^j_u^ and thus one gets a 'renormahzation' of the 
spinning particle's mass: 

/i = — . 30 

Remark 2. The first space-time curvature of the particle's world line governed 

by (22) is constant, i.e. -;— ll'itll = 0. 

Remark 3. The point symmetries of (20) are being exhausted by pseudo-ortho- 
gonal (resp. conformal) transformations if m 7^ (resp. m = 0). The proof may be 
found in [10]. 

We can present two different (a = 1, 2) Lagrange functions which produce the 
equation (20), 

^(«)=(l + v.v)V.(l+g..||vAe(.)P)^ -Ml + v.v) . (31) 

These differ by the total time derivative: 

d 

L(2) ~ ^(i) = -J- arctan ■ 



dt VI + VaV^ 

Now it's time to pass over to the hamiltonian counterpart exposition. 

3. HAMILTONIAN DESCRIPTION 

A detailed exposition of the (generalized) hamiltonian theory applicable to the 
odd-order differential equations at no less extent than to the even-order ones may 
be found in [6]; a concise exposition is presented in [7]. 

A conjectural Legendre transformation is given by the momenta 



dL 


d dL 


d\i 


dt dM' 


dL 




dM' 





The Hamilton function 

H = -L + p.\i + p'.v' 

in case of a regular Legendre transformation gives rise to the following Hamilton 

system of first order equations: 

dH dy'd d 

"Pfe-TlPa = 0, [62) 



dx" dx^dt" dt 

dH d „ , ^ 

OPa dt 

dH dv^ d . (dy" dv^\ d , 

+ ^3lx'+hz7-Tr7 hlP. = 0. (34) 



ap'„ dxt'dt V-^Pfe dp'^) dt 
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Both Li and Lj from (31) produce one and the same Hamihon function: 

*(v Av') 



H 



;i + vv)3/2 



^(1 



v-v 



-1/2 



= p.v + /iVl + V-V. (35) 

This happens because both L^ and Lj lead to the same 'zero-order' momentum, 

*v' V 



'l+Vv)3/2 



fl 



VT 



v-v 



(36) 



However, the 'first-order' momentum p' has the identically equal to zero first (resp. 
second) component if one starts with L^ (resp. Lj) alone. This suggests that we take 

in an attempt to proceed further with a kind of regular Legendre transformation. In 
this case one calculates out the following expressions for the two components of the 
momentum p': 



2Vl + vv(l+vivi) 



P2 



2Vl + V-V(l + V2V2 



(37) 



In (34) the inverse to the Jacobi matrix of the Legendre transformation is in- 

dispensable to know (at least in its jr-i part). The Jacobi matrix itself is readily 

obtained from (36) and (37), so we cite below only that part of its inverse, important 
for the forthcoming calculations. 



dy"" 



I V2V^(3 + 3V2v2-F2ViV^) 



|(^ 



v-v 



,3/2- 



X 



-1 



'1 + V2V 



2\2 



ViV^(3 -|- 3 ViV^ + 2 V2V^ 



V 



'1 +ViV 



1\2 



(3J 



/ 



where A denotes the determinant of the matrix in (38). The Reader may also easily 
convince himself by direct calculation that 



<9p b 



0. 



(39) 



As far as the Hamilton function (35) does not depend on the space and time 
variables, the essential part of the Hamilton system, the equation (32), constitutes 
nothing else but the conservation of the momentum. 



^ = 0. 

dt 



(40) 
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The remaining two equations mean merely that we pick up only the holonomic 
solutions of (40). Thus, the equation (33) in view of (39) reads 

dx 

and so it describes the kernel of the first one of the two differential forms in (8). The 
equation (34) in view of (38) reads 



dp' A \ dt 

a TT 

But we can calculate -^— ^ from (35) by means of (38) and also differentiate the 

definitions (37) explicitly. After substituting into (41) it becomes evident that the 
Hamilton equation (34) describes the kernel of the second differential form in (8). 
In fact the equation (41) transforms into the following system: 

M . (rfv - v'dt) = , 

where M is the matrix from (38). 

Let us return to the gravitational dipole particle. The Dixon system of equations 
in the general curved four- dimensional space-time accepts arbitrary parametrization 
of the particle's world line. It reads: 



d^ 
ItX 
d& 



d (42) 

2^Au. (43) 



The force ^ depends, apart from space and time variables, also upon the particle's 
tensorial spin & and velocity u; it must satisfy the constraint ^.u = o . Once 
supplemented by the Pirani constraint (4), the system (42, 43) allows an equivalent 
transcription in terms of the spin four- vector (29), 'resolved' with respect to the 
variable ^. Namely, in place of (43) we can write (renouncing the constraint ||u|| = 
1) 

*P = ^u- (sgng) * (iiAuAs), (44) 

S Au = O. 

In the parametrization by time (u" = 1), denoting by P the three- vector part of the 
four-vector ^, the definition (44) takes up the shape 

P = mo^=^ + — ^f^(v'xs-5»v'xv), (45) 

Vl + v-v (l + v-v)-v^ 
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where we use the notation v for the three-vector part of the velocity four-vector u 
in the special case when u" = 1; also s = {s°, s) according to the general template. 
Reassuming the particle's motion be two-dimensional only, one obtains from (45) 

It suffices to compare (46) with the expression (36), assuming the notation (30), to 
come immediately to the conclusion that the quantity P coincides with 5^ times the 
canonical hamiltonian momentum. 

Proposition 3 The Hamilton function for the planar motion of free relativistic 
spinning particle with the mass (30) is given by (35). 

Remark 4. The Mathisson equation in general-relativistic framework and under 
arbitrary world line parametrization was considered in [11]. A set of the Lagrange 
functions to produce it in fiat space-time was suggested in [12]. A second-order 
differential-geometric connection related to the equation (20) was constructed in [13]. 
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